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Abstract 

X-rays diffraction pictures time dependence of deformed alloy Pd-Ta being charged with hydrogen 
has been shown to be possibly caused by multi-pits character of energetic relief in the states space. 
Phenomenological model representing alloy structural evolution as a occasional roaming on minima 
of internal energy of non-ergodic system, has been offered in Lorenz synergetic scheme frames. Here, 
order parameters are the part of minima occupied by the system, conjugated field is considered 
to be Edwards entropy and control parameter is taken to be internal energy. Thermodynamics 
interpretation of Pd-Ta-H alloy evolution structure as a complex non-ergodic system is offered 
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1 Introduction 

The main object of statistics physics has lately become complex non-ergodic systems - spin and structural 
glasses, non-ordered heteropolymers, granulated matter, transportation flow, etc (see Q,). The main 
peculiarity of such systems is phase space separated into isolated regions with every region representing 
meta-stable thermodynamic state and their number A'o exponentially exceed full number of quasi particles 
N ^ oo, N: No — exp (sN) , where s is so-called Edwards entropy respective to every complexity |H, 
Unlike Bolsman's measure describing disorder in given statistics ensemble, Edwards entropy describes 
disorder in complex system states distribution on internal energy minima respective to statistics ensemble. 
Statistics ensemble's play the role of particles in complicated systems. The distribution on these ensembles 
is described with effective temperature T and entropy S, introduced by Edwards So, complicated 
system with T = is corresponding to granulated matter distributed eventually on all minima (flat 
distribution) . 

Edwards systems that have exponentially huge number of minima A'^ ^ 1, give us an example which 
is opposite to ergodic systems where Nq = 1. According to this, an intermediate case where the number 
of minima Nq is neither exponentially big nor equal to one is of great interest. 

Deformed Pd-Ta alloys being electrolytically charged with hydrogen represents this example as it is 
shown below. Really, solid state material under the influence of strong deformation can go so far away from 
the equilibrium state that the self-organization effects acting in dissipative structure formation become 
essential ^. During hydrogen charging of metal its atoms create essential internal tensions that can 
lead the system into a non-equilibrium state as an external influence does. The researches of annealed 
Pd-W-H alloy § deformed Pd-Sm-H, Pd-Er-H, Pd-Mo-H and Pd-Ta-H alloys § |, |, |l^, |l], |l|, g 
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show that the non-monotonous structure transformation of non-regular character can take place there. 
Using of synergetic models |l^ allows to explain the main features of these systems behavior. But 
these models give no explanation of peculiarities of Pd-Ta-H alloy structural evolution shown below (see 
also [p2l). This paper is devoted to this problem solution. 



2 Experimental data 

The method of the investigations was described in details in |l^, |l^, |2| , so we point out only at the 
fact that the original charging of Pd-Ta deformed sample {7at.%Ta) has been done electrolytically for 
15 minutes at current density of 40 mA/cm2. After relaxing for 176 hours, the sample has been charged 
for the second time for the same period of time at 80 viA/m^ . X-Rays researches of diffraction maxima 
220 311 were done with computerized diffractometer with using CuKai components of x-rays spectrum 
doublet. The diffraction maxima decomposition into components was completed with "Origin" program 
package taking into consideration Lorenz form of components ||lj, [l^ . 

The led researches show that, after the hydrogen charging of Pd-Ta deformed alloy, regular displace- 
ment of diffraction maxima which represents the lattice getting bigger and after that smaller take place. 
During the relaxation, stochastic changes of diffraction maxima position, their broadness and symmetry 
appear. So does stochastic transformation of their form including appearance and disappearance of some 
picks. Diffraction curves and their decompositions received for 200 maxima after the secondary charging 
at the different moments of time are shown at the fig.l. They represent complex structure of diffraction 
picks and non-monotonous character of their changes with time. Actually, the diffraction maxima have 
the structure shown at the table here. 

As a bell-like form of diffraction maxima reflects homogeneous phase, one can presume that such 
changes of diffraction pictures are consequences of phases mutual transformations accordingly to Pd, Ta 
and H atoms re-distribution over different regions of the system. These structural changes can be repeat- 
ing but not periodic. On the other hand, different broadness of components and their non-monotonous 
time changes mean that this process is accompanied with non- monotonous defect structure evolution 
and elastic fields created by this evolution appearance. Along with essential difference between hydrogen 
binding energy with ideal and defect Palladium lattice M], the existence of far-acting fields causes es- 
sentially non-equilibrium character of researched system As it is known from the theory of glass- like 
theory |jl[, it causes complicated energetic relief appearance in the phase space of non-equilibrium state. 
Character of this relief is defined by the original defect structure and alloy atoms distribution, that are 
formed by mechanical work up, regimes of hydrogen charging and following relaxation. The difference of 
time dependence for the 220 and 331 coherent scattering regions oriented differently to the surface of the 
sample points out to energy being passed from freedom degrees to others that can lead to diffusion flows 
turbulence |p^ . 



3 Theoretical scheme 

The non-monotonous structural changes discussed above can be explained with the preposition that Pd- 
Ta deformed alloy hydrogen charging leads to the set of meta-stable states appearance corresponding 
to different phases and defect structures. In this case, generalized system diffusion over internal energy 
minima in accordance with the pointed states leads to stochastic structural changes observed during the 
experiment. 

To understand this phenomenon, one needs to take into consideration the non-monotonous character 
of evolution, which is not of periodic character. These changes remind of strange attractor behavior 
described with synergetic Lorenz scheme [|l^. Earlier, using Lorenz scheme has allowed us to describe 
non- monotonous behavior of Pd-Er-H double- phased alloy M. While doing that, we parameterized the 
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system with value portion of phase enriched with Er, concentration of Er atoms there and Er atoms traps 
concentration in matrix. 

This parameterization cannot be used in muhi-phased Pd-Ta-H alloy where the situation is much 
more difficult. Facing this problem for the first time, using parameters generalized over all the phases or 
yielding the most typical one helps avoiding the necessity of many phases behavior description. But this 
way we could have lost the extremely important peculiarity of multi-phased Pd-Ta-H alloy evolution: not 
only value correlation of the phases but their number change with time. According to this, the problem 
of parameters choice for the synergetic model which tries to describe Pd-Ta-H system stochastic behavior 
becomes very important. 

According with Ruelle-Takens theorem [Q, non-trivial behavior of such system is observed if the 
number of parameterizing freedom degrees is over two. At phase transition, the system behavior is 
described with hydrodynamic mode and its amplitude represents order parameter defined with thermostat 
state. Self-organized (synergetic) system peculiarity is: not only thermostat's influence over the subsystem 
but also its influence over the thermostat is very important. This influence can be either direct or indirect. 
First of them is defined with conjugated field, the second one - with control parameter. So, for the alloys 
of ordered synergetic model , order parameter is in accordance with usual long range order parameter, 
conjugated field refers to difference between chemical potential of components and control parameter - to 
difference of unit cell knots population of different atoms Usual dissipative regime of phase transition 
is realized if relaxation time for order parameter and conjugated field is much longer then it is for control 
parameter. Otherwise, stationery value of control parameter increasing over crucial value leads not to 
ordering of this system but to transforming into regime of strange attractor ]20| . 

To explain non-monotonous change of Pd-Ta-H alloy we take into consideration the fact that hydro- 
gen charging leads to internal energy complicated relief state formation in the space. This relief has 
many minima separated with barriers (multi- valley structure). This preposition usage explains critical 



slowing down of structural transformation in the Pd-H systems [21 . This slowing down is caused by the 
hierarchical portrait of the relief; big minima are overlapped with smaller ones, then the smaller minima 
are covered with the minima that are even smaller, etc. Because of that, during its evolution, system is 
enforced to fill the smallest minima first, starting to fill bigger and bigger minima up to the biggest one 
which describes the system as a whole subject. 

As to the observed case, the consideration which is limited with a type of minima that are not 
essentially different from each other is enough for our purpose. We numerate them with index a covering 
the set of numbers 1,2,.. Nq. As a result, the system evolution is described with probability of distribution 
over these minima Pa changes. Defining of Pa (t) time dependence requires Fokker-Planck equation 
to be solved which is very difficult |^^ . However, for the preliminary stage, we consider the integral 
characteristics that describe the type of Pa distribution and thermodynamical behavior of the system as 
a whole subject enough for investigation. We take parameters as following: 

• Number N of internal energy minima where the system is situated at the moment (obviously, it 
describes a half of dispersion width for probability Pa ) ; 

• Entropy of the system dispersion over those minima 

S = -^Palnpa (1) 

a 

• Specific internal energy at given minimum a, which defines the full value of the internal energy 

'paCa (2) 



We will use the relation n = N/Nq of number N of internal energy minima occupied by system and 
full number of minima in the system Nq as an order parameter. Then specific Edwards entropy s = S/Nq 
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which defines disorder in the distribution over minima a is described with conjugated field and internal 
energy density and acts as a control parameter. 

For the phenomenological description of the system evolution, main parameters ^ of velocity 

changes are to relate to their values n, s, e. Here, Lorenz system advantage is in its accordance to the 
simplest choice of Hamiltonian in the respective microscopic representation |^ . Linear Lorenz equation 
is 

dn n 

-77 = \-9nS, (3) 

at Tn 

where first part describes Debuy's relaxation of n(t) value to n = with character time t„; second 
part which contains positive coefficient g„ describes number of internal energy minima occupied by the 
system at increasing of entropy s in the dispersion over those minima. 

Unlike equations describing velocities ^ contain non-linear parts that represent the influence 
of the subsystem over the thermostat, showing the feedback we have discussed before. Thus, the equation 
for entropy change is 

ds s , ^ , 

-T7 = \-gsne, (4) 

dt Ts 

where first part describes Debuy's relaxation to s = with time r^. In accordance to second principle of 
thermodynamics, non-linear part which containes positive coefficient takes into consideration entropy 
s increasing because of system dispersion over less deep minima of internal energy. 

Last of the required equations defines velocity of internal energy increasing. 

de e - eo , ... 
— = \-gesn (5) 

at Te 

Here, Debay's relaxation with character time Te leading not to zero value of the internal energy e 
but to Co defined by system position at the phase diagram and the alloy preliminary work up is taken 
into consideration. Non-linear part with positive coefficient represents negative feedback meaning 
that the re-distribution over internal energy minima should lead to decreasing of its full value e because 
of the system transformation to deeper minima. It should be pointed out to a very important role of 
non-linear parts. Their competition describes self-organized system behavior: positive feedback in 
causes entropy s increasing because of relation between relative number of minima n and internal energy 
e. Positive feedback in equation (^) describes internal energy e decreasing because of the correlation 
between the number of minima n and entropy s. 

The system of differential equation (^)-(|^) generally has no analytic solution. Using of non-dimensional 
values n, s, e, t scaling them as shown below seems to be very suitable: 

n,n = (rsgs)^ (rege)^ , s,„ = em = {T„gn)~^ {Tsgs)~^ , T„. (6) 



As a result, the equations (||)-(||) appear to be 



- = + (7) 



ds 

^dt + 
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de 

= (i? - e) - sn, (9) 
where the non-dimensional parameters are entered 

T = ^, e = — , E^^. (10) 

In adiabatic regime r, <C 1 the principle of hierarchic co-ordering takes place. In accordance with 
this principle, changes of conjugated field s(t) and control parameter e(t) follow order parameter n{t) 
changes In this case, left parts of the equations (|^), (^) can be omitted and they lead to: 

(U) 



In accordance with thermodynamics principles, they mean that Edwards' entropy monotonously grows 
but internal energy falls down with phases number growing up. On the other hand, variable n excluding 
from the equalities (Bl]) gives simple dependence of entropy on internal energy: 



s = y/e (E-e). (12) 

Using the definition of the temperature 



E\ ^ [e 



we get 

It means that Edwards' temperature becomes negative in the interval of the internal energy values 
E/2 < e < E and the system is in the process of self-organization. Therefore, in accordance with 
definition (^3|) at T < every increasing of internal energy Se > leads to entropy decreasing Ss < 0, 
i.e. to self-ordering. 

Substitution of the first of the equalities (^ into equation (^ turns Lorenz' system into Landau- 
Khalatnikov equation. 

dn dW , , 

where the role of the free energy is played by synergetic potential 



W =^n^ -^lii{l + n'^) (16) 

At little values of internal energy E which is saved by the system as a result of external influence, 
dependence W(n) monotonously increases with minima at n = and multi-phased state is not advan- 
tageous synergetically. However, exceeding over crucial value E = 1 synergetic potential W(n) reaches 
minima at 

no = (17) 
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Therefore, dissipative process of the system relaxation into the state respective to the value ng 7^ of 
relative number of the phases appears to be advantageous. During this process, Edwards' temperature 



(113) reaches stationary value 



T„ = -^^, (18) 



and , being negative at over-critical values E > 1, monotonously decreases with growing of the internal 
energy saved by the system in consequence of external influence. 

The represented adiabatic regime r, <C 1 reflects monotonous transition of the dissipative system 
into non-ergodic stationary state. As we are interested in non-monotonous behavior, we can make a 
conclusion that the character time r, Q ratio should not be small. If it happens, analytical investigation 
of Lorenz' system appears to be impossible but the analysis |po| show that the regime of strange attractor 
of the self-organizing system specific for the experimental situation is realized at the condition 8 > r > 1. 
Thus, time of internal energy change should exceed respective values of entropy and number of minima 
occupied by the system. Obviously, it is realized in the experiment investigated Jl^ . 

Besides these limitations for the ratio of relaxation time, strange attractor existence requires execution 
of i? > 1 condition |^ . This means that the stochastic changes of the structure in consequences of system 
states re-distribution over the minima respectively to different phases are realized only at the condition 
that the value of the internal energy saved by the system in a result of external influence exceeds critical 
value Cm which is deflned with the third expression of (^). It can be noticed, that the systems that have 
great value of positive feedback's linear constant gn and non-linear constant gs and long relaxation time 
T„, Ts are predisposed to stochastic behavior. 



4 Conclusion 

The led analyses shows that the stochastic structure transformation is realized because of the competitive 
influence of order parameter over the conjugated held and control parameter of self-organizing system. 
In the case of multi-phased Pd-Ta-H alloy where the roles of two last parameters are played by specific 
value of entropy s and internal energy e this competition is supported with the basic thermodynamic 
equality: 

f = e-Ts (19) 

where T - Edwards' effective temperature. This equality is a result of Edwards' temperature definition 
( |l3| ) and conjugated entropy definition. 

.-IL.. (20, 

One of us has shown recently |Q that the pointed definitions along with identity ( p^ ) - are originally 
from the simplest field scheme in the limitation of which self-organizing system behavior is parameterized 
with double-component fields of order parameter, conjugated field and control parameter. Here first 
components are reduced to the values used before-respective number of occupied minima n, entropy s 
and internal energy e of non-ergodic system. Second components represent generalized flows conjugated 
to the pointed values such as flow of the probability density q re-distribution among minima of internal 
energy; thermodynamics force —V/ equal to gradient of free energy, taken with the opposite sign; tem- 
perature opposite gradient — VT. These components change appears to be essential, if the system is at 
non-stationary state. Its example represents non-monotonous behavior of multi-phased Pd-Ta-H alloy. 



6 



Hopefully, the synergetic scheme offered here represents the variant of strongly non-equilibrium state 
systems thermodynamics theory and its development is far from being completed [ p4| , p5[ . 

According to the analysis shown above, kinematics condition of the system of non-monotonous behav- 
ior consists of internal energy changes and temperature gradient changes that exceed respective value of 
entropy and gradient of free energy as well as of minima occupied by the system and conjugated probabil- 
ity flow. This condition explains abnormally big values of time intervals where stochastic behavior appear. 
Dynamic condition requires entropy increasing to have essential influence over the increasing number of 
minima where the system is situated . On the other hand, it requires decreasing of internal energy as 
consequences of re-distribution in order to lead to essential increasing of the entropy itself. Obviously, 
these conditions are provided with the essential values of gn and parameters that are realized at not 
high barriers separating minima and low "rigidity" of the internal energy relief. Thus, thermodynamics 
condition is considered a requirement of strongly non-equilibrium state of the system. This state of the 
system is reached with preliminary work up of the system such as radiation, essential deformation, hard- 
ening, etc. For the Pd-Ta-H system being investigated this work up consisted of preliminary deformation 
and hydrogen charging that lead to elastic tensions that are about 10% of elastic module character value. 
Considering this reasoning, we can presume that the internal energy saved by the system is big, too. 

Let's point out at two observations of general character. First one is- critical slowing down of the 
structural transformation found in the system Pd-H can be also found in the Pd-Ta-H alloy considered 
here. Our preliminary data received after hydrogenation for several times confirms this presumption. 
This means that at the quantity of hydrogenations essential increasing, the considered above multi- valley 
structure which is represented with the same minima transforms into multi-leveled hierarchical structure 
of energetic relief, which leads to the system evolution slowing down. This slowing down was found by 
the authors of Q. 

Second observation is about the data of multi-phased deformed Pd-Mo-H alloy where the non- 
monotonous structural changes that are analog to the ones considered above were found |Q. A s for 
Pd-Ta-H regular displacement of diffraction maxima takes place. This process reflects the lattice parame- 
ters getting bigger and, after that, smaller in anisotropic way. During the following relaxation, stochastic 
changes of the diffraction maxima components number, their positions, their width and intensity occur. 
Using this data, one can conclude that the non-regular behavior of the strongly dissipative multi-phased 
system is a specific feature for the structural changes of the alloys after hydrogen charging. 
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Time after the charging, h 
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Number of figure 
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Table 1: Time dependence of diffraction maxima (220) number 
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Figure 1 : Time change of the position and profile of the diffraction line (220) after second charging: a) 1 
h; b) 3,5; c) 6; d) 11; e) 25,5; f) 28; g) 30,5; h) 33; i) 73,5; j) 97,5; k) 148 (vertical line shows the position 
respective to original condition - 173 h of relaxation after second charging) 
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